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Proton  transport  resistance  in  the  polymer  electrolyte  of  proton  exchange  membrane  fuel  cells  (PEMFCs) 
is  a  significant  contributor  to  ohmic  overpotential.  Conductivity  is  largely  controlled  by  water  uptake, 
which  requires  a  delicate  water  balance  to  avoid  an  increase  in  mass  transport  overpotential  due  to  dif¬ 
fusion  resistance  in  porous  components.  The  present  work  explores  this  relationship  with  distributed 
liquid  water  and  high  frequency  resistance  (HFR)  measurements.  The  HFR  response  to  local  changes  in 
liquid  water  content  is  characterized  for  two  gas  diffusion  layers  (GDLs)  with  different  thermal  resis¬ 
tivities.  From  these  experiments,  it  is  shown  that  liquid  water  must  be  present  in  the  gas  diffusion 
layer  (GDL)  to  minimize  the  HFR.  Furthermore,  the  impact  of  varied  saturation  levels  induced  by  GDL 
thermal  gradients  on  drying  transients  is  considered  with  both  in-plane  and  through-plane  neutron 
imaging. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

PEMFC  water  management  is  a  broad  research  topic  because 
of  the  wide  range  of  length  scales  present  in  a  fuel  cell  stack. 
Researchers  must  consider  transport  from  the  molecular  scale 
in  the  polymer  electrolyte,  to  droplet  formation  in  the  dif¬ 
fusion  layer  pore  structure,  to  small  slugs  of  water  in  the 
bipolar  plate  channels,  to  larger  slugs  (on  the  scale  of  centime¬ 
ters)  accumulating  in  the  bipolar  plate  manifolds  and  ancillary 
equipment.  If  water  was  not  required  to  facilitate  proton  con¬ 
duction  between  the  anode  and  cathode  catalyst  sites,  operating 
conditions  could  easily  be  prescribed  to  enable  removal  of 
all  product  water  in  the  vapor  phase.  However,  the  current 
state-of-the-art  in  membrane  materials  clearly  indicates  that  first- 
generation  PEMFC  systems  with  acceptably  high  performance 
at  reasonable  cost  will  require  a  humid  operating  environ¬ 
ment  that  results  in  water  condensation  at  low  stoichiometric 
ratios. 

The  relationship  between  liquid  water  accumulation  and  trans¬ 
port  resistance  has  been  widely  studied  with  both  direct  and 
indirect  diagnostic  methods.  In  a  recent  review  of  electrochemi¬ 
cal  and  physical  diagnostic  methods,  Wu  et  al.  presented  many  of 
these  techniques  and  the  conclusions  based  on  them  [1,2].  From  a 
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commercial  stack  perspective,  the  majority  of  water  transport 
research  is  related  to  the  water  activity  required  for  membrane 
hydration.  The  relationship  between  the  proton  transport  resis¬ 
tance  and  operating  parameters  has  been  widely  studied  with 
impedance  in  the  high  frequency  spectrum  [3-5].  Liquid  water, 
which  can  block  gas  transport  (observed  via  direct  visualization 
techniques),  can  be  problematic  during  operation  at  high  relative 
humidity  (RH)  -  a  condition  which  is  required  to  minimize  pro¬ 
ton  resistance  [6].  The  trade-off  between  too  much  water  (high 
oxygen  transport  resistance)  and  too  little  water  (high  proton  trans¬ 
port  resistance)  results  in  optimum  performance  near  a  water 
balance  independent  combined  anode  and  cathode  exhaust  stream 
RH  of  1 00%,  as  maximum  membrane  conductivity  is  achieved  while 
avoiding  severe  mass  transport  loss  associated  with  excessive  liq¬ 
uid  water  [7]. 

Combining  visualization  with  local  current  and  impedance  mea¬ 
surements  has  enabled  researchers  to  investigate  the  relationship 
between  liquid  water  accumulation  and  local  performance  loss 
[8-10].  These  experiments  universally  demonstrate  that  severe 
water  accumulation  negatively  impacts  local  current  output.  Many 
researchers  generalize  any  liquid  water  content  in  the  carbon  fiber 
substrate  as  ‘flooding’,  which  has  a  negative  connotation.  However, 
these  studies  did  not  consider  the  optimum  saturation  (the  fraction 
of  pore  volume  in  the  compressed  substrate  that  contains  liquid 
water),  as  the  presence  of  some  liquid  water  might  be  necessary  to 
reduce  proton  transport  resistance  to  a  minimum  value  for  a  given 
membrane  and  ionomer  type.  If  this  is  the  case,  the  gas  diffusion 
layer  (GDL)  should  be  designed  to  maintain  a  low,  but  non-zero, 
level  of  saturation. 
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Fig.  1.  Distributed  measurement  tool,  (a)  printed  circuit  board,  (b)  composite  flow 
field  plate  bonded  to  board. 


In  the  current  work  we  use  the  thermal  resistivity  of  the  GDL 
to  impact  water  saturation  level  and  correlate  these  changes  to 
the  proton  transport  resistance  in  the  membrane.  Expanding  on 
previous  work  [8,9],  we  have  developed  an  unique  combination  of 
experiments  and  analyses  to  consider  the  proton  transport  resis¬ 
tance  associated  with  liquid  water  accumulation,  as  it  is  known 
that  proton  resistance  in  perfluorosulfonic  acid  solid  electrolyte  is 
only  minimized  if  some  liquid  water  is  present  [11,12].  By  moni¬ 
toring  distributed  high-frequency  resistance  (HFR)  and  local  liquid 
water  content  during  transient  drying,  this  article  considers  the 
saturation  level  required  to  maximize  ionic  conductivity  in  the 
membrane,  and  how  it  relates  to  the  magnitude  of  saturation  in 
the  GDLs. 

2.  Experimental 

The  experimental  work  conducted  in  this  study  combined  a 
small-scale  fuel  cell  with  distributed  HFR  measurements  and  neu¬ 
tron  imaging.  A  fuel  cell  test  stand  with  individually  calibrated 
control  parameters  was  used  to  operate  test  cells.  The  PEMFC  test 
fixture  was  specifically  designed  within  commercial  stack  con¬ 
straints  to  provide  pure  counter-flow  in  the  active  area  where 
anode  and  cathode  gas  flow  is  parallel  in  opposite  directions.  This 
design  is  described  in  detail  by  Owejan  et  al.  [9]. 

2.  t.  Current  distribution  measurement 

The  fuel  cell  current  distribution  tool  consisted  of  a  printed 
circuit  board  with  isolated  contacts  that  redirect  current  through 
shunt  resistors.  The  circuits  were  printed  on  a  polytetrafluoroethy- 
lene  (PTFE)  substrate  to  improve  neutron  transmission  in  a  similar 
manner  to  work  reported  by  Gagliardo  et  al.  [8].  As  shown  in 
Fig.  la,  the  measurement  array  consisted  of  25  columns  down-the- 
channel  with  4  rows  across  the  channels  and  lands.  The  voltage 
response  across  each  shunt  was  measured  with  a  separate  data 
acquisition  channel.  The  flow  field  was  machined  in  thin  graphite 
composite  resin  and  bonded  to  the  distributed  measurement  board 
(Fig.  lb).  The  continuous  flow  field  results  in  slight  “blurring”  of  cur¬ 
rent  between  segments  but  this  will  be  no  larger  than  the  1.5  mm 
thickness  of  the  plate.  The  plate  was  verified  to  be  isopotential  by 


Table  1 

Measured  diffusion  and  thermal  parameters  for  the  GDL  material  studied  [13-15]. 


Parameter 

TGP-H-060 

Grafil  U-105 

Unit 

Thermal  resistivity  [13] 

3.1  ±0.3 

9.2  ±0.3 

cm2  KW_1 

Effective  diffusion  coefficient  [14] 

0.036  ±0.002 

0.037  ±0.002 

cm2  s-1 

Uncompressed  thickness 

225  ±12 

230±  10 

fxm 

Electrical  Resistivity  [15] 

7.8 

13.3 

m£2  cm 

Porosity 

66% 

73% 

- 

monitoring  voltage  along  the  perimeter  at  varying  distances  from 
an  applied  current. 

To  measure  the  high  frequency  impedance  component  associ¬ 
ated  with  the  zero-phase  resistance,  the  current  was  perturbated 
at  a  frequency  of  1  kHz  with  a  peak-to-peak  amplitude  of  500  mA. 
At  this  high  frequency,  the  voltage  and  current  were  assumed  to  be 
in  phase  and  ohmic  resistance  was  calculated  based  on  the  voltage 
response,  assumed  constant  for  all  segments  (isopotential  plate) 
relative  to  the  individual  current  response  which  varied  based  on 
the  resistance  of  each  segment. 

2.2.  Materials 

Two  GDL  substrates  with  different  thermal  conductivities  but 
similar  diffusion  resistances  were  selected  to  demonstrate  how 
through-plane  thermal  gradients  impact  liquid  water  content  in 
the  GDL.  By  varying  the  GDL  thermal  conductivity  the  fuel  cell 
water  content  could  be  altered  at  a  given  operating  condition, 
thus  establishing  the  correlation  between  liquid  water  quantity  and 
membrane  resistance.  The  high  and  low  thermal  conductivity  GDLs 
selected  for  this  study  were  Toray  TGP-H-060  and  Mitsubishi  Rayon 
Corp.  (MRC)  Grafil  U-105,  respectively.  Each  sample  was  prepared 
with  7  wt%  PTFE  in  the  substrate  and  a  30  [xm  thick  microporous 
layer  consisting  of  an  8:1:1  carbon: PTFE :FEP  ratio.  The  relevant 
measured  parameters  for  the  resulting  GDLs  are  listed  in  Table  1. 
The  GDLs  were  compressed  to  20%  strain  based  on  the  measured 
uncompressed  thickness.  All  experiments  were  run  with  an  18  |xm 
thick  membrane  and  0.2/0.3mg  Ptcnrr2  anode/cathode  catalyst 
loading  membrane  electrode  assembly  (MEA)  manufactured  by 
W.L.  Gore. 

3.  Results  and  discussion 

The  GDLs  being  compared  in  this  study  were  expected  to 
exhibit  different  water  accumulation  behavior  at  higher  temper¬ 
ature  operating  conditions  due  to  the  temperature  gradient  that 
exists  between  the  coolant  fluid  and  the  cathode  catalyst  layer 
[16].  As  water  condenses  in  the  porous  medium,  the  temperature 
gradient  will  dictate  the  location  of  the  condensing  liquid  front. 
Hence,  a  lower  temperature  gradient  will  result  in  more  liquid 
water  volume,  all  else  held  constant  (diffusion  properties,  for  exam¬ 
ple).  This  behavior  was  confirmed  with  in  situ  neutron  imaging 
at  the  National  Institute  of  Standards  and  Technology  (NIST)  Neu¬ 
tron  Imaging  Facility  (NIF)  by  measuring  the  integrated  liquid  water 
thickness  through  all  layers  with  the  neutron  beam  perpendicular 
to  the  active  area  plane  [9].  Shown  in  Figs.  2  and  3  is  a  comparison  of 
liquid  water  content  and  performance  for  the  Toray  and  MRC  GDLs 
at  various  current  densities  during  a  polarization  curve  acquired  at 
80  °C  cell  temperature.  The  similar  polarization  performance  with 
a  nearly  factor  of  two  difference  in  liquid  water  volume  indicates 
that  GDL  water  content  can  have  an  insignificant  influence  on  fuel 
cell  performance  for  certain  operating  conditions. 

The  data  shown  in  Fig.  4  at  80  °C  demonstrate  the  contrast  in  liq¬ 
uid  water  accumulation  due  to  the  nonlinear  relationship  between 
saturation  vapor  pressure  and  temperature.  As  stated  above,  even 
though  there  is  a  factor  of  approximately  two  difference  in  water 
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Fig.  2.  Comparison  of  liquid  water  accumulation  in  MRC  and  Toray  GDL(80°C,  1. 5/2.0  H2/air  stoichiometry,  95%  inlet  RH,  200  kPa  abs.  outlet  pressure,  counter-flow). 
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Fig.  3.  Performance  comparison  of  Toray  and  MRC  GDLs  (80  °C,  1. 5/2.0  H2/air  stoi¬ 
chiometry,  95%  inlet  RH,  1501<Pa  abs.  outlet  pressure,  counter-flow). 


content  at  all  current  densities  for  80  °C  cell  temperature,  the  volt¬ 
age  and  HFR  measurements  are  nearly  identical  (Fig.  3).  Only  at  the 
1.5  A  cm-2  condition  is  there  some  evidence  of  greater  mass  trans¬ 
port  loss  in  the  cell  with  high  thermal  conductivity  GDL  material. 
To  further  demonstrate  that  water  accumulation  in  the  bulk  of  the 
GDL  is  mainly  related  to  the  induced  temperature  gradient  between 
the  cathode  catalyst  layer  and  coolant,  the  liquid  water  content 


at  35  and  80  °C  is  presented  in  Fig.  4.  For  each  of  these  measured 
water  distributions,  the  associated  temperature  gradients  were  cal¬ 
culated  based  on  the  thermal  resistivities  in  Table  1  and  the  heat 
flux  at  1.5  A  cm-2  determined  from  the  overpotentials  in  Fig.  3.  At 
35  °C  a  difference  of  1.5  kPa  in  saturation  vapor  pressure  gradient 
results  from  temperature  gradients  of  2  and  6  °C  for  Toray  and  MRC 
papers  respectively.  Conversely,  at  80  °C  the  difference  in  induced 
temperature  gradient  for  the  two  materials  results  in  a  difference 
of  8.8  kPa  in  saturation  vapor  pressure  gradient.  Consequently,  the 
total  water  content  for  the  Toray  material  is  less  dependent  on  cell 
temperature.  Flowever,  the  in-plane  liquid  water  content  along  the 
channel  direction  is  clearly  influenced  in  Fig.  4.  At  80  °C,  the  liquid 
water  is  more  concentrated  at  the  center  of  the  active  area  while  the 
anode  and  cathode  inlet  regions  remain  relatively  dry.  This  temper¬ 
ature  dependence  is  later  used  to  consider  variations  in  GDL  water 
content  at  the  same  operating  condition. 

This  relationship  between  accumulated  water  and  performance 
was  explored  further  by  monitoring  local  FIFR  relative  to  liquid 
water  content  for  the  cell  assembled  with  MRC  GDL  (Fig.  5).  Flere, 
the  neutron  image  was  analyzed  over  discrete  areas  aligned  with 
the  segmentation  of  the  current  distribution  tool  shown  in  Fig.  1. 
Each  water  thickness  measurement  is  the  integrated  liquid  water 
measured  per  pixel  that  is  averaged  over  an  area  equivalent  to  the 
current  measurement  segment.  The  plot  shown  in  Fig.  5  demon¬ 
strates  the  strong  dependence  between  bulk  water  in  the  GDL  and 
membrane  proton  resistance  that  is  governed  by  its  hydration  level. 
For  this  plot,  the  gradient  in  water  thickness  values  was  achieved 
by  preconditioning  the  cell  at  a  wet  condition  (35  °C,  0.4  A  cm-2, 
1 50  kPa  outlet  pressure,  2/2  H2/air  stoichiometric  ratios)  then  purg¬ 
ing  at  a  drier  exit  condition  (35  °C,  0.1  A  cm-2,  101  kPa  outlet 
pressure,  2/16  H2/air  stoichiometric  ratios).  At  this  low  temper¬ 
ature  condition,  the  cathode  purge  dries  the  cell  slowly  (~10min), 
thus  allowing  several  levels  of  liquid  water  thickness  to  be  captured 
with  30  s  neutron  integration  times  in  each  segment.  The  inset 


Fig.  4.  Impact  of  GDL  induced  thermal  gradient  on  measured  water  volume  at  35  and  80  °C  (0.4  A  cm-2, 1 .5/2.0  H2/air  stoichiometry,  95%  inlet  RH,  1 50  kPa  abs.  outlet  pressure, 
counter-flow). 
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Fig.  5.  HFR  response  to  local  water  thickness  during  a  slow  drying  transient;  liquid 
equilibrated  equivalent  liquid  water  thickness  in  the  ionomer  is  indicated  by  the 
vertical  dashed  line.  Neutron  images,  binned  water  thickness  and  HFR  distributions 
that  represent  data  are  also  shown.  Precondition:  35 °C,  0.4 A  cm-2,  150 kPa  abs. 
outlet  pressure,  2/2  H2/air  stoichiometric  ratios;  drying  condition:  35  °C,  0.1  A  cm-2, 
101  kPa  abs.  outlet  pressure,  2/16  H2/air  stoichiometric  ratios. 
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Fig.  6.  Drying  transient  comparing  HFR  response  for  Toray  and  MRC  GDLs.  Pre¬ 
condition:  80°C,  0.4Acirr2, 150 kPa  abs.  outlet  pressure,  2/2  H2/air  stoichiometric 
ratios;  drying  condition:  80°C,  0.1  Acm-2, 101  kPa  abs.  outlet  pressure,  2/16  H2/air 
stoichiometric  ratios. 


neutron  images  and  distributed  HFR  in  Fig.  5  demonstrate  this  dry¬ 
ing,  as  a  front,  moving  from  the  cathode  inlet  to  outlet  during  the 
length  of  the  purge.  These  experiments  were  run  under  a  wide  array 
of  conditions,  and  the  same  general  trend  in  HFR  was  observed  as 
integrated  liquid  water  thickness  decreased. 

In  Fig.  5,  the  sharp  increase  in  HFR  as  water  thickness  approaches 
a  value  that  is  associated  with  liquid  equilibrated  ionomer  water 
content  (an  estimated  13  pm  equivalent  water  thickness  for  the 
950  equivalent  weight  membrane  at  20  moles  water  per  acid  site) 
is  somewhat  obvious  given  the  strong  dependence  of  the  ionomer 
“lambda”  value  (number  of  water  molecules  per  sulfonate  site)  on 
the  phase  of  water.  However,  what  is  unique  about  these  data  is  that 
the  HFR  asymptote  indicates  only  a  small  amount  of  liquid  water 
in  the  GDL  is  required  to  maintain  the  ionomer  lambda  value  at  a 
liquid  equilibrated  state. 

This  difference  in  water  storage  capacity  realized  by  GDL 
thermal  properties  can  have  varying  impacts  on  fuel  cell  sys¬ 
tem  performance.  On  one  hand,  a  higher  stored  volume  of 
liquid  water  can  buffer  high  load  and  temperature  transients. 
On  the  other  hand,  a  lower  stored  volume  can  shorten  the 
freeze  preparation  purge,  thereby  reducing  hydrogen  consump¬ 
tion.  We  explore  this  further  by  comparing  the  local  HFR  response 
of  the  two  GDLs  selected  for  this  study  during  a  drying  tran¬ 
sient  at  80  °C.  Data  are  presented  for  segments  of  the  HFR 
distribution  board  at  both  upstream  and  downstream  locations, 
1.4  and  16.9  cm,  respectively,  from  the  cathode  inlet  edge  of 
the  active  area.  The  HFR  response  in  Fig.  6  clearly  demon¬ 
strates  this  trade-off  where  the  upstream  and  downstream  HFR 
transients  are  much  slower  for  Toray  relative  to  MRC  GDLs. 
With  regard  to  purge  energy  this  delays  the  rapid  increase  in 
HFR  (indicative  of  substantial  drying  in  the  MEA)  for  the  less 
thermally  resistive  Toray  GDL  by  as  much  as  25  s  at  the  down¬ 
stream  location,  which  translates  to  a  significant  efficiency  loss 
as  compared  to  the  MRC  GDL.  In  practice,  the  longer  shut¬ 
down  purge  duration  required  for  the  Toray  material  could 


be  balanced  by  the  benefit  of  increased  robustness  to  long 
excursions  at  dry  operating  conditions.  The  sensitive  relationship 
shown  in  Fig.  6  is  important  when  selecting  GDL  materials  for  a 
specific  application  or  to  optimize  a  specific  region  of  the  operating 
space. 

Finally,  the  integrated  water  measurement  presented  thus  far 
is  ignoring  the  anode  vs.  cathode  distribution  in  these  two  experi¬ 
mental  configurations.  To  elucidate  this  we  applied  high  resolution 
neutron  imaging  methods  to  measure  the  through-plane  water 
content  by  orienting  the  cell  such  that  the  neutron  path  was  parallel 
to  the  active  area  plane  [17].  This  comparison  is  provided  in  Fig.  7, 


Through-Plane  Thickness 


Fig.  7.  High  resolution  neutron  imaging  of  through-plane  water  profiles  for  Toray 
(80  °C,  0.6  A  cm-2,  4/4  H2/air  stoichiometry,  95%  inlet  RH,  150kPa  abs.  outlet  pres¬ 
sure,  counter-flow). 
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where  integrated  liquid  water  thickness  is  scaled  to  GDL  saturation. 
Here  it  is  observed  that  the  increased  overall  water  content  with 
the  Toray  material  is  a  result  of  both  a  change  in  the  anode  distribu¬ 
tion,  where  more  water  is  stored  toward  the  anode  flow  field  and 
an  increase  in  cathode  GDL  liquid  water  content  with  a  similar  dis¬ 
tribution.  This  result  indicates  a  shift  in  the  overall  water  balance, 
which  correlates  to  the  slower  HFR  response  of  Toray  observed  in 
Fig.  6,  as  the  additional  anode  water  removal  by  the  cathode  gas 
stream  is  slowed  by  the  membrane.  Also,  an  increase  in  anode  water 
further  complicates  the  trade-offs  mentioned  previously,  because 
draining  water  from  the  anode  subsystem  can  impact  system  effi¬ 
ciency. 

4.  Conclusions 

The  relationship  between  GDL  saturation  and  proton  trans¬ 
port  resistance  was  examined  with  local  measurements  of  liquid 
water  and  HFR  by  considering  GDLs  with  different  thermal  con¬ 
ductivities.  Experiments  comparing  Toray  with  MRC  GDLs  yield 
similar  polarization  performance  to  1 .5  A  cm-2  while  the  3  x  higher 
thermal  resistance  MRC  material  retains  significantly  less  liq¬ 
uid  water  at  80  °C.  This  difference  in  water  content  at  the  same 
operating  condition  was  used  to  consider  how  local  GDL  sat¬ 
uration  impacts  proton  transport  resistance  in  the  membrane. 
Results  show  that  this  relationship  is  very  sensitive,  and  a 
small  amount  of  liquid  water  in  the  GDL  is  required  to  min¬ 
imize  proton  transport  resistance,  but  as  saturation  increases 
beyond  this  value,  the  proton  transport  resistance  remains  con¬ 
stant.  Through-plane  liquid  water  distribution  and  transient  drying 
experiments  highlight  the  practical  trade-off  associated  with 
this  relationship.  A  lower  thermal  resistance  GDL  can  result  in 
higher  anode  and  cathode  GDL  liquid  water  content,  provid¬ 
ing  a  water  buffer  for  dry  transients  but  also  increases  required 
purge  energy  for  reliable  cold  starts  in  automotive  applica¬ 
tions. 
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